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The statistical mechanics of a linear non-interacting polymer chain with a large number of 
monomers is considered with fixed angular momentum. The radius of gyration for a linear polymer 
is derived exactly by functional integration. This result is then compared to simulations done with a 
large number of non-interacting rigid links at fixed angular momentum. The simulation agrees with 
the theory up to finite size corrections. The simulations are also used to investigate the anisotropic 
nature of a spinning polymer. We find universal scaling of the polymer size along the direction of 
' the angular momentum, as a function of rescaled angular momentum. 
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Polymer chains in a vacuum have been shown to have great importance in many experimental situations, most 
I ■ often in application to mass spectroscopy techniques that utilize the desorption of proteins into a vacuum. The ability 
I to desorb very large molecules without compromising their integrity has many uses in biology [l| . Understanding the 
statistical mechanics of such systems may help to improve current spectroscopic techniques by providing information 
1 1 > ' about other quantities, aside from mass, that may be usefully measured. Aside from other applications 0, there is 
C "^j also the intrinsic interest in understanding such systems. 

The dynamics of vacuum polymers has been the subject of a number of recent investigations 043 an d show many 
unusual features different from what are seen in solution. Polymers close to the coil-globule transition show small 
Lyapunov exponents Q. However in many situations, the dynamics are in accord with detailed microcanonical calcu- 
lations, implying that these systems are ergodic @. Ideal chains show highly oscillatory time correlation functions Q, 
whereas the addition of self interactions damps these oscillations although they are still quite pronounced for small 
chains }5jj . Although there is no coupling to a heat-bath, the nonlinear dynamics of these models gives rise to an effec- 
tive damping of individual monomers which is similar to that of Kelvin damping Q . Adding local angular potentials 
to polymers increases this type of damping but does not inhibit oscillations for sufficiently long chains [J] . 

In this paper, we further investigate the equilibrium statistics of a noninteracting polymer chain of N links in a 
vacuum, where energy, momentum, and angular momentum are conserved. Previous work Q found an analytical 
solution for a ring polymer with conserved angular momentum. The average radius of gyration as a function of angular 
momentum was obtained. The most surprising feature of that calculation is that the radius of gyration when the 
angular momentum is zero is less than that of a chain without that conservation law. Although the radius of gyration 
. still is proportional to y/~N but now with a different proportionality constant. 

The calculation performed for a ring chain Q does not extend easily to linear chains. Here we use an approach 
00 ■ that diverges considerably from the previous one to obtain the statistics of a linear chain under the same conditions. 
In both cases, although the intermediate steps produce rather lengthy expressions, the final answers are quite simple. 
This suggests that there may be another underlying principle that can be used to understand systems with constant 
angular momentum. 

Aside from the average radius of gyration of such a system, there is another important quantity that describes 
it, that has hitherto not been studied. For finite angular momentum, the chain is expected to become anisotropic, 
flattening in the direction of the angular momentum vector. It is therefore of interest to calculate the average size of 
a chain both perpendicular and parallel to the angular momentum. However our analytic method is only applicable 
in calculating quantities which have rotational symmetry such as the radius of gyration. It will not work for studying 
chain anisotropy. Therefore we turn to numerical methods to study this problem. We find, rather surprisingly, that the 
chain continuously flattens in the direction of the angular momentum showing universal scaling in the same rescaled 
variables used to characterize the full radius of gyration. 

II. EXACT SOLUTION 

Following previous work Q we begin by writing down the entropy as a function of the total energy E, angular 
momentum L, and number of particles N, including terms for the conservation laws to be enforced. 
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where K is the total kinetic energy and $ the total potential energy. We are choosing a coordinate system so that 
the center of mass r cm = 0. 

We take the potential energy to be the elastic energy of the polymer chain and in addition include an isotropic 
quadratic potential that will be used to calculate the radius of gyration. In terms of the continuous arclength s 
dependent position variables that we will use below, 

< 5 *=OS |2+er2< * ),ia ' <2> 

where I is the step length. 

Because we are taking N to be large, we can convert this to a partition function, with the energy conserving delta 
function expanded in exponential form. 
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After integrating over the momenta and taking into account rotational symmetry we arrive at the form Q 
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Z[P,L,N] = y sin(kL)C(l3,k)dk (4) 
L Jo 

where £ is a function of the magnitude of k and inverse temperature (3 only. The constant c is of no physical importance 
and 

C(P,k) = /e-r , ^ +e H-W !, +vW a )+«W !, +*l*l"«to, (l . c|n)5r(a)> (5) 



Now we perform a scaling on L, k, and s in order to normalize out the length of the chain and the temperature. 



L< = — — ^=L (6) 

k' = i^^k (7) 
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After doing these rescalings, the form of Eq. 2] remains unchanged, and Eq. [5] changes to 



(8) 



C(ftfc') = / e"X <*"/M+.')h'<.') I +»'<.')")+.'='(-') 2 +*l^l , ^ a(l .; m)Jr ' (s . ) (9) 



It is interesting to note that the explicit dependence on j3 has disappeared and ( is now only a function of k! . The 
dependence on the total length of the chain has also been factored into k', but the length of the individual monomers 
remains. Now we expand the center of mass conserving delta function in complex exponential form. 



C(fc')= / e-/o^ (fc ' 2/2 ' +e ' )( *' 2+y ' 2)+eV2 +^^ (10) 



These terms can be combined with the existing exponents using completion of squares which results in path inde- 
pendent phase factors and a constant offset to the paths, e.g. x'(s) — > x'(s) + const. . Since all possible paths are 
included in the integral, the offset has no effect. 

We now have a nice separation into a simple Gaussian integral and a factor that is very similar to the path integral 
form of a harmonic oscillator. 



The first integrals of the Z's are trivially integrated but we must still perform the path integral in parentheses. These 
decouple into path integrals over the x, y, and z directions and each one has an action of the form 



/ — (x 2 +wlx 2 )dt. (12) 
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M does not depend on k or e, and u z has no k dependence. After making this analogy, we can see that for any two 
fixed path endpoints, the functional integration will give the Greens function for a thermal harmonic oscillator. To 
obtain a sum over all paths, we simply need to integrate the Greens function p [§| over both its arguments. 

Ct = J Pi(x 1 ,x 2 )dxidx 2 (13) 

where i = 1,2,3 represent the x, y, and z, path integrals respectively. The Greens function is easily integrable and 
the result is 

Ci « 1 (14) 

J uji cosh (-\/l2wj) 

Putting all of this back into Eq. [4] we get the form for Z 



Z(L') = -J %= / fc'sin(fc'L') ^^dk'. (15) 

L\j sinh( v / T2w z ) Jq sinh(v / 12w) 



The radius of gyration can be expressed as 



2 _ 1 dlnZ 

R v--M— U =° (16) 



The logarithmic derivative allows us to separate the factors in Z and treat each as terms in a sum. The first term is 
simple 

e^o de V sinh (yi2w z ) 18 V 7 

The second term is more complex due to the integral, so we begin by rewriting the logarithmic derivative as a ratio. 

^Ck'sm(k'L') . " dk' 

lim , / Slnh( ^ ) (18) 

e->0 ^k'smtk'U) -ft- dk' V 7 

Since the integral is over k' we can pull the derivative and limit inside the integration. This simplifies the integrand 
considerably as lim e _>.o^ = fc'/v3 

# [ n °°/c'sin(fc'L')— -JH—-dk' r> N 2 
-0 J^k'Mk'L^ ^^ dk' 6 ^3tanh(^) 

Put all of the pieces together to obtain our final result. 

M 2 18 67rtanh(^) V ' 

A plot of this is shown in Fig. Q] by the solid line. The simulation results also displayed will be discussed in Sec. 
El Note that for V = 0, R 2 g /{Nl 2 ) = (1/18) + 2/(3vr 2 ) w 0.1231, which is about 74% of 1/6, the value of the same 
quantity without conservation of anguar momentum enforced. 

For the same reason as for the ring calculation [6], for fixed L' and N — > 00, we are permitted to take this system 
to be at constant temperature, because as was discussed in detail [6J, the relation between energy and temperature 
has no dependence on L' in this limit. 



III. HIGH L LIMIT 



In the limit of large L, Eq. I2TJ1 gives R 2 /Nl 2 — > L'/6ir. One would expect typical configurations at high L would 
consist of a fairly straight line rotating rapidly around its center. The solution for the equilibrium case of constant 
angular velocity can be found by setting the change in tension along the chain equal to the centripetal force. This 
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FIG. 1: Plot of -j^z along with simulation results using a chain length N — 128. Inset is a section from L' — 5.5 to L' = 7 
showing the exact solution as a solid line along with results from simulations done with N — 32 (farthest from line) and N = 64. 



gives an equation for the position of each monomer as a function of its location along the chain, r(s), for simplicity 
defined with s = at the center of the chain. For large N, the tension approaches zero at the ends of the chain. 



(21) 



where k is the entropic elastic spring coefficient. From this we get that r(s) oc sin(y / ^ws), and that Nloj/2^/k — tt/2. 
Higher modes exist, but they double back and are expected to not be minimal free energy solutions. They are also 
not the straight chain model we are assuming. Using the integral definition of angular momentum 
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the radius of gyration squared is 
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And we see that in the limit of high i, our model behaves as one would expect. 
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FIG. 2: Plot of P(L )), with clear non-Gaussian behavior for high L' 



IV. PARTITION FUNCTION 

The partition function itself can be obtained easily from Eq. [HI by taking e — > and performing the integration. 

«-»o LJ sinh(2fc') 32L' C osh 2 (L'f) V ' 

To obtain the distribution over L' we need to normalize this Z, taking in to account the angular integration over the 
direction of L' , 

rllrrjTl 2 



so our normalized partition function is 



^L' A Zdl! = n z (25) 



z 1 tanh ( £, f) f26) 



32L' cosh 2 (i'f 

A plot of this is displayed in Fig. [2] The distribution shows a long exponential tail into high V 



V. SIMULATION RESULTS 



In addition to the exact calculation, further studies were carried out using simulation. First the analytical results 
concerning the radius of gyration as a function of Li were verified. Second we were able to investigate other quantities 
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FIG. 3: Plot of radius of gyration in the direction of angular momentum as a function of rescaled L. 



that were beyond the means of our analytic method. 

We used a molecular dynamics method [J] that was developed to simulate chains in a vacuum. It consisted of freely 
rotating rigid links, conserving energy and angular momentum. Rigid links were employed to minimize problems with 
equilibration that are often seen with one dimensional nonlinear systems [ill Il2| . The input angular momentum 
and the output measurements were scaled by the total chain length for comparison. The initial values of the angular 
momentum were chosen to cover a range of L' values, and the angular momentum was explicitly checked and conserved 
during the runs. First the simulation was compared to the theoretical results for the normalized radius of gyration. 
The results of these simulations can be seen in Fig. [TJ and show excellent agreement with the theoretical prediction. 
The main plot shows the exact result (solid line) along with data for chains with N — 128. The inset shows that 
the deviation in the asymptotic form for high L' decreases as the number of simulation chain links is increased and 
appears to be due to the finite size of the simulation system. The two chains lengths used in the inset are N = 32 
and TV = 64 which are more accurate than the data for N = 128. They also indicate that the finite size corrections 
to the analytic form are 0(1/N). 

For nonzero angular momentum the polymer chain is expected to become anisotropic. The analytic methods used 
earlier require an isotropic form for the quantities being averaged, causing the investigation of the anisotropy by such 
means to be outside the scope of our analysis. Therefore we turned to the simulation to explore the aspect ratio of 

the polymer as a function of L' The aspect ratio is defined as «/ 1 (z 2 ) / (Rf) .The aspect ratio itself is dominated by the 

asymptotic linear behavior of the radius of gyration, and displays what appears to be an inverse power law falloff as 
expected. More interestingly, (z 2 ) itself appears to fall off as L' increases. This is shown in Fig. [3] where the vertical 
axis is scaled in the same manner as for the radius of gyration, by dividing by and taking I, m, T all to be unity, 

and the horizontal axis is L' The three chain lengths seem to collapse onto one curve after the rescaling of length and 
angular momentum. This is unexpected because for a ideal gaussian chain each dimension has independent statistics. 
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In the limit of high L the rigid link model should approach a straight line solution. However this is unlikely to be 
an explanation, as the straight line regime would imply a leveling off of the radius of gyration, which is not seen. 
Moreover, the falloff is independent of the number of chain links in the simulation and collapsed onto a single curve, 
so is not likely due to the non gaussian nature of the model. This decrease is slow and does not appear to have a 
nonzero asymptote, and is likely to be a low power law or logarithmic in nature. 

VI. CONCLUSIONS 

In this paper, we considered the equilibrium properties of an ideal linear chain in a vacuum that conserves energy, 
momentum, and angular momentum. We were able to compute the average radius of gyration of such a chain as a 
function of its angular momentum L. We also computed the distribution of angular momenta for chains in thermal 
equilibrium. We verified that our analytical result for the radius of gyration is correct by performing numerical 
simulation and by analyzing its asymptotic form in the limit of large angular momentum. The derivation of this 
result differs from that of a ring chain but in both cases the final result is relatively simple involving hyperbolic 
trigonometric functions. The underlying reason for this is still unclear. 

Our numerical simulations show that the radius of gyration perpendicular to the angular momentum vector increases 
with L, as to be expected, however in addition to this, the radius of gyration parallel to the angular momentum, 
which we take to be in the z direction, decreases. This is very different than what a naive analysis would suggest. If 
we were to go to a frame rotating at angular velocity 51, then in thermal equilibrium (without angular momentum 
conservation), this is equivalent to an additional potential — mO rSL/2, where r xy is the projection of a coordinate 
onto the x — y plane. For an ideal gaussian chain, all three directions decouple and(z 2 ) is independent of L. However 
a more rigorous analysis based on the approach of this paper is not equivalent to this, and it is not clear how statistics 
in the x — y plane and z direction are coupled to cause this effect. 
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